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Abstract 

We present spin wave dispersions for MnO, NiO, and a-MnAs based on the recently-developed 
quasiparticle self-consistent GW method (QSGM^), which determines an optimum quasiparticle 
picture. For MnO and NiO, QSGVF results are in rather good agreement with experiments, in 
contrast to the LDA and LDA-I-C/ descriptions. For a-MnAs, we find a collinear ferromagnetic 
ground state in QSGTV, while this phase is unstable in the LDA. 



PACS numbers: 71.15-m,71.10-w,71.20-Eh 



I. INTRODUCTION 



The magnetic linear response is a fundamental property of solids. It is given by the 
spin susceptibility when the spin-orbit coupling is neglected (as we will do in this article). 
The spin susceptibility is equivalent with the spin fluctuations, as can be seen from the 
fluctuation-dissipation theorem. Low-energy spin fluctuations can control some low energy 
phenomena such as magnetic phase transitions, and contribute to resistivity through spin-flip 
scattering of electrons. Antiferro(AF)-magnetic spin-fluctuations can play an important role 
in high-Tc superconductoraii^. It is also the central quantity entering into the description of 
quantum-critical phenomenai^i^. We expect that reliable first-principle methods to calculate 
the spin susceptibility should give important clues to understand these phenomena. In this 
paper, we concentrate on the magnetically-ordered systems, where the spin susceptibility 
should be dominated by spin waves (SW) at low energy. 

In spite of the recent development of such methods, we still have large class of systems 
where we can hardly calculate the spin . A typical 

example is MnO; Solovyev and Terakura gave an analysis for the calculation of its SW 
energies^. Then they showed main problem is in the non-interacting one-body Hamiltonian 

from which we calculate the non-interacting spin susceptibility used for the calculation of 
the SW energies. given by the local density approximation (LDA), LDA+f/, or even the 
optimized effective potential (OEP)^ are not adequate. In the LDA+U case, they traced the 
error to a misalignment of the 0{2p) bands relative to the Mn(3c/) bands. It is impossible 
to choose the U parameter to correct the misalignment, because the U parameter can only 
control the exchange splitting within 3d bands. A possibility may be adding some other 
parameter in addition to U so as to correct the misalignment; however, such a procedure 
including more parameters become less universal. This situation is somehow similar with the 
case of optical response (dielectric function) calculation for semiconductors, where given 
by LDA is with too small band gap, thus requires some additional correction like scissors 
operator. Our case for the spin susceptibility for MnO is rather worse; LDA supplies too 
problematic to be corrected in a simple manner. 

Another possibility is to obtain by some hybrid functional; it has been shown that it 
can work as explained below, however, it could be problematic from the view of universality. 
Muscat, Wander, and Harrison claimed that a functional BSLYP-^^ (containing 20% of Fock 
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exchange) works even for solids. However, Franchini, Bayer, Podloucky, Paier and Kresse^ 
showed that a httle different functional, PBEO is better than B3LYP in order to obtain better 
agreement with experiments as for the exchange interaction. PBEO is a combination of 25% 
of the Fock exchange with a generalized-gradient approximation(GGA)^i. However, such a 
functional could be not so universal, mainly because the effect of screening (therefore the 
ratio of the Fock exchange) are dependent on materials. In fact, de P. R. Moreira, Illas, and 
Martini^ reported that a hybrid functional containing 35% of Fock exchange gives best results 
for NiO; the ratio of the Fock exchange is rather different from the case of MnO by Franchini 
et al. This is somehow consistent with latest careful examinations by Fuchs FurthmuUer, 
Bechstedt, Shishkin and Kresse^^ and Paier, Marsman, and Kresse^; they clarify the fact 
that a hybrid functional should be limited, because the screening effects (corresponding 
to the ratio of the Fock exchange) can be material-dependent. These seem to indicate a 
difficulty to pick up an universally-applicable hybrid functional. This difficulty becomes 
more problematic when we treat inhomogeneous systems, e.g, to treat the Schotkey-barrier 
problem, where the screening effects are very different in metal side and in semiconductor 
side. 

Considering these facts, it is necessary to start from good without such problems. Our 
recently-developed quasiparticle self-consistent GW method (QSGW^), which includes the 
above screened exchange effects in a satisfactory manner— li^ii^ii^ii^i^SiSi. QSGiy determines 
a reference system of representing optimum quasiparticle (QP) picture in the sense of 



19, it is based on a self-consistent 



Landau-Silin Fermi liquid theory. As discussed in Reft 
perturbation theory within all-electron full-potential GW approximation {GWA), but it is 
conceptually very different from the usual full self-consistent GW . QSGVT self-consistently 
determines not only , but also the screened interaction W ^ and the Green's function G 
simultaneously. 

We have shown that QSGiy gives QP energies, spin moments, dielectric functions and so 
on in good agreement with experiments for wide range of materials. There are systematic 
but a little disagreement from experiments. For example, as shown in Fig.l in Ref ItI, we 
see error that calculated band gaps are systematically larger than those by experiments. 
A recent development by Shishkin, Marsman and Kresse^i confirmed our conjecture^'' that 
the inclusion of electron-hole correlation effect in W will correct the error. Their method 
is a simplified version of the full Bethe-Salpeter equation (BSE) for W] it includes only the 
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static and spatially-local part of the first-order term in the BSE, based on the procedure 
given by Sottile, Olevano, and Reining^. However, since the band-gap error itself is small 
enough, such simplifications may cause no problem. Considering this success, we believe 
that QSGW is a basis for future development of the electronic structure calculations. 

In this paper, we treat MnO and NiO with AF ordering II (AF-II)2^, and a-MnAs. 
a-MnAs is NiAs-type grown on GaAs epitaxially, thus is a candidate for spintronics 
applications^i. For this purpose, we have developed a procedure to calculate the spin sus- 
ceptibility at zero temperature. It is a general procedure for a given self-consistent method 
which determines H^, even when contains non-local potentials as in the Hartree-Fock 
method. We then apply it to LDA, and to QSGW. After we explain the method in the 
next section, we will show SW energies obtained with QSGW are in good agreement with 
experiments for MnO and NiO. See RefJlOl for dielectric functions for NiO and MnO. For 
MnAs, our calculation shows that a collinear FM ground state is stable in QSGW though 
it does not in LDA. 

At the end of introduction, we give a discussion to justify using the one-particle picture 
(band picture) of "Mott insulator" for MnO and NiO; it is essentially given by Terakura, 
Williams, Oguchi, and Kiibler in 1984^'^ (in the following discussion, "charge transfer type" 
or "Mott type" does not matter). Based on the one-particle picture, the existence of some 
spin moment (or exchange splitting, equivalently) at each cation site is very essential to 
make the system insulator. This is consistent with the experimental facts that all the 
established "Mott insulator" are accompanied with the AF (or some) magnetic ordering. 
Thus the concept "Mott insulator vs Band insulator" often referred to is misleading, or 
rather confusing. In order to keep the system insulating, any ordering of spin moment is 
possible provided the system retains a sufficiently large enough exchange splitting at each 
site (we need to use the non-coUinear mean- field method). In this picture, metal-insulator 
transition at zero-temperature (e.g, consider a case to compress NiO) is nothing but the 
first-order transition from magnetic-phase to the non-magnetic phase described by a band 
picture. On the other hand, the transition at finite temperature to para-magnetic phase 
occurs because of the entropy effects due to the accumulation of SWs; then the transition 
is not accompanied with the metal-insulator transition because the exchange splitting (or 
local moment) at each site is kept even above the Neel temperature T^. This picture is 
very different from that assumed in Refsl5,25, where they emphasize the priority of their 
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method LDA + U + "dynamical mean field theory (DMFT)". On the contrary to their 
claim, we insists that our treatment should be prior and much closer to reality for such 
systems, because of the following reasons. 

i) One-particle treatment in our QSGW allows us to perform parameter-free 
accurate calculations where we treat all the electrons on the same footing; this 
is very critical because of the relative position of cation 3d bands to 0{2p) 
is important (also their hybridization; we have no SW dispersion without hy- 
bridization). Further, we are free from uncontrollable double-counting prob- 
lemM, nor the parameter like U which is externally introduced by hand. In con- 
trast, LDA+f/+DMFT carries these same problems which are in LDA+[/, or 
rather highly tangled. Thus it is better to take a calculation by LDA+f/ +DMFT 
as a model in cases. As an example, we guess that the distribution probability 
of the number of 5/ electrons in 5-Plutonium calculated by LDA+[/+DMFT— 
will be easily changed if we shift the relative position (and hybridization) of 5/ 
band with respect to other bands. 

ii) The DMFT at zero temperature takes into account the quantum-mechanical 
onsite fluctuation which is not included within the one-particle picture; it allows 
a system to be an insulator without magnetic order. However, we expect that 
such quantum-mechanical fluctuation is not essentially important to determine 
its ground state for materials like NiO and MnO. This is based on our findings 
that QSGW results can well reproduce the optical response^^i^, and also the 
magnetic responses as shown in this paper. These QSGW results are not perfect, 
however, supphes us a good enough starting point. For example, in order to 
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describe the d-d multiplet intra transitions (e.g, see Fig.6 of Refj! 
and Minami; they are very weak in comparison with interband transitions), it 
may be easier to start from the cluster models or so; however, parameters used 
in these models will be determined by QSGWeven in such a case. 

iii) At finite temperature, the DMFT can take into account not only such 
quantum-mechanical fluctuations, but also the onsite thermal fluctuations si- 
multaneously; this is an advantage of DMFT. However, in MnO and NiO, low- 
energy primary fluctuations are limited to the transverse spin fluctuations except 
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phonons. These can be included in DMFT but it is essentially described by the 
local-moment-disorder— as the thermal average of the one-particle picture. Thus 
no advantage of DMFT if only the thermal fluctuations are important. 

II. METHOD FOR SPIN SUSCEPTIBILITY CALCULATION 

We may divide first-principle methods to calculate SW energies into three classes; (A), 
(B), and (C). (A) is from the Heisenberg Hamiltonian, whose exchange parameters J are 
determined from the total energy differences of a set of different spin configurations^^i^ii^^. 

(B) and (C) are based on perturbation. (B) estimates J from static infinitesimal spin rota- 
tionsS^i^i. We go through the Heisenberg model even in (B). In contrast, (C) determines SW 
energies directly from the poles in the transverse spin susceptibility (r, r', t — t') (defined 
below) in the random phase approximation (RPA) or time-dependent LDA (TDLDA)^'^'^. 

(C) gives the spectrum including life time, and spin-flip excitations. Because (C) is techni- 
cally difficult, (A) or (B) have been mainly used. (B) is regarded as a simplification of (C); 
but real implementations entail further approximations. 

Our method belongs to (C). Our formalism is applicable to any even if it contains 
non-local potential. At the beginning, we introduce some notations to treat the time-ordered 
transverse spin susceptibility 

X+-(r,r',t-t') = -^(T(5+(r,t)5-(r',tO))- (1) 

(...) denotes the expectation value for the ground state; T(...) means time-ordering, and 
S^{r,t) = S^{r,t) ±iS'^{r,t) are the Heisenberg operators of spin density. Since we as- 
sume coUinear magnetic ordering for the ground state, we have {S^{T,t)) = {S^{r,t)) = 0; 
2{S^{r,t)) = M(r) = n\r) — n^{r). n\r) and n^{r) mean up and down electron densi- 
ties. Ma{r) is the component of M(r) on the the magnetic sites a in unit cell. The Fourier 
transform of is 

X+-(T + r,r» = ^Ee^"^xr(r>r',^), (2) 

q 

where T is a lattice translation vector, and the number of sites, r, r' are limited to a unit 
cell. 
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Next we derive two conditions Eq. (jlj) and Eq. ([5]) below, which rigorously satisfies. 
Taking the time derivative of Eq. ([1]), we obtain 



d_ 
di 



- 1 rfVx+-(r',r,t'-t) 

d'r'{T{[[H, 5+(r', t')], ^-(r, t)])) - z / rfV([5+(r', t'), S'^r, t)])6{t' - t), (3) 



where [A, B] = AB — BA. H denotes the total Hamiltonian of the system. We have used 
^^"^dt' ^ ~ i[H , S^{r\t'y\. We assume H has rotational symmetry in spin space, so that 
[ff, / d?r'S^{r' , t')] = 0. Then the first term in the right-hand side is zero. The second term 
reduces to M(r) because [S^{r',t), ^"(r, t)] = 2S'^(r,t)5(r — r'). Thus Eq. ([3]) is reduced to 
be 

/^rfVx+-o(r',r,u;) = ^, (4) 

where fl denotes the unit-cell volume. Note that Eq. (jlj) is satisfied for any u. At u 0, 
this means that M(r) is the eigenfunction of Xqro(r, r', cj) with divergent eigenvalue; this 
is because a magnetic ground state is degenerate for homogeneous spin rotation. Another 
condition is the asymptotic behavior as — oo. It is given as 

X+-(r', r, u) ^ ^5(r - r') + 0{l/u;'). (5) 

00 

This can be easily derived from the spectrum representation of We use Eq. and 

Eq. (j5]) to determine the effective interaction U in the following. 



As in Ref. 



38l . we define the effective interaction U (r, r', uj) as the difference between 



(x+ ) -"^ and the non-interacting counterpart: ) ^ {r,r',uj) 



-1 



(x+-)-^ = U"+- +U. (6) 



In TDLDA, U is the second derivative of the exchange- correlation energy, f/(r, r') = 
—6'^Exc/SS^{r)6S^{r') = /xc(r)(5(r — r'), which is local f/(r,r') oc 6{r — r'), c<j-independent, 
and positive. Then we can show that in TDLDA satisfies conditions Eq. (jl]) and Eq. ([5]) 
automatically^-. In the case of containing nonlocal potentials (e.g. in the case of the 
Hartree-Fock method), U is no longer independent of u. This is because the natural ex- 
pansion of x~^~ in the many-body perturbation theory requires solving the Bethe-Salpeter 
Eq. for the two-body propagator r2ir3, r4, cj), thus U defined in Eq. is not di- 



rectly identified as a kinds of diagrams. Ref, 



371 did not pay attention to this point. We can 



calculate x' 



.0+- 



in Eq. as 



,0+ 




(7) 



where k' = q + k. (r, r',ci; = 0) is negative definite matrix. Our definition of and 
also can be different in sign from other definitions in the literature because we start 

from Eq. ([1]). 

In order to realize an efficient computational method, we assume that the magnetization 
is confined to magnetic atomic sites, and we explicitly treat only a degree of freedom of spin 
rotation per each site. Then we can determine U with the help of Eq. (jlj) and Eq. ([5]) as in 
the following. As a choice to extract the degrees of freedom, we consider a matrix D{q,uj) 
as 



and i5*^(q, C(j) defined in the same manner. The dimension of the matrix D{c[,uj) is the 
number of magnetic sites. Here we define ea(r) = Ma{r)/Ma where Ma = Jad^rMa{r); and 
define ea(r) so that ea(r) oc ea(r) and J d^rea{r)ea{r) = 1; thus ea(r) = ea(r)/ (iV(ea(r))^. 
Corresponding to Eq. IQ, we define the effective interaction {U{q,uj))aa' as 



For the calculation of /^"(q, uj) from Eq. ([7]), we use the tetrahedron technique^^, which allow 
us to use fewer k points in the first Brillouin zone (BZ) than those required for the sampling 
methocP^. U defined in Eq. ([9]) should include all the downfolded contributions from all 
the other degrees of freedom. We now assume that U is q-independent and site-diagonal, 
so that it can be written as f/aa/ (q, cj) = Ua{uj)6aa'- Since Eq. (jl]) reduces to a constraint 
J2a'{D{q = 0,ci;))a'a = Ma/uj, we determine Ifaiuj) from 






(9) 




With this Ua{uj) for Eq. (Q, we finally have 



(D(q,^)) 



1 



(10) 



a 




(11) 
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Eq. (ISl) reduces to (-D(q, uj))~,l -^^aa' at ^ oo; (-D(q, uj))~}^ given by Eq. ([10]) gives this 
correct asymptotic behavior. Note that we determine U just from the requirement Eq. (jl]) 
because of our approximations "onsite only [/" and "a basis per magnetic site" . If we need to 
go beyond such approximation (e.g. muhiple basis per site), it will be necessary to introduce 
additional informations, e.g. a part of x{^.i = 0) evaluated by numerical linear-response 
calculations (perform the QSGVT self-consistent calculations with bias fields). By Fourier 
transformation, we can transform (D(q, Ci;))a'a into DfiRi{uj)] the same is also for D^, J and 
so on. Here R = Ta is the composite index to specify an atom in the crystal. For later 
discussion we define 

J(q,.;) = -(D°(q,^))-^ + -^, (12) 

where D^^^u) is shorthand for -DxaXal'^)) it is T independent. The second term in Eq. ( IT2|) 
is included just in order to remove the onsite term from J. Then Eq. f|TT]) can be written as 

J(q,w) = J(q,^) - I — jSaa'- (13) 

Here, the second term (onsite term) in Eq. ( fT2l) is irrelevant because of the cancellation 
between two terms in Eq. (fT3ll . 

The preceding development for (Z}(q, c<j))^^ facilitates a comparison with the Heisenberg 
model, whose Hamiltonian is H = — J2r J2r' Jrr'^R'^r' = Ta). As shown in Appendix 
lAt the inverse of the susceptibility in the Heisenberg model is: 

P^(q,^))-^ = ^5aa'- J^(q), (14) 

where Ma = \2Sji\. Let us compare Eq. (fT^ with Eq. (fTOl) . This J^/(q) is given by 
Eq. flAgp . which is almost the same as Eq. f|T3|) : only the difference is whether we use or 
J. This suggests how to construct the Heisenberg model which reproduces Eq. ( fTOl) as good 
as possible; a possibility is that we simply assign J(q, u = 0) (neglecting the cj-dependence) 
as J^(q). We have confirmed that this approximation is good enough to reproduce SW 
energies in the case for MnO and NiO. However, it is not true in the case of a-MnAs; 
then we have used another procedure given by Katsnelson and Lichtenstein^^: we identify 
J{c[,uj = (SW energy at q)) as J^. This construction exactly reproduces SW energies 
calculated from D'^''^. 



9 



As a further approximation to calculate J(q, c<j = 0), we can expand it in real space as 
(omit uj for simplicity) 

R R R R' 

where we use Eq. f lT2|) : we use notation that onsite part Z)^ = -D^jj and the off-site part 
^r!r' — ^%R' ~ ^RR'^RR'- Here we have used the assumption that -D^^/ are small in com- 
parison with onsite term D^. This approximation corresponds to the usual second-order 
perturbation scheme of the total energy; if the spin rotation perfectly follows the rotation 
of the one-particle potential, 7^- is trivial; it is equal to the difference of the one-particle 
potential between spins (exchange-correlation potential in the case of DF) because is 
the inverse linear-response to determine the one-particle potential for given spin rotation. 



331134. In cases, this approxi- 



Essentially the same equation as Eq. ( |T5l) was used in Refs. 
mation is somehow mixed up with the "long wave approximation" to expanding J around 
D°(q = 0)^^; however, they should be differentiated. In order to have rough estimate of Jrr', 
we can further reduce this to the two sites model as originally presented by Anderson and 
Hasegawa^i^. For an AF magnetic pair (half- filled case), we obtain the following estimate: 

Jnu'-^—^D,,,— ^-^^, (16) 

where t denotes the transfer integral, and AE^^ is the onsite exchange splitting. We have 
Dr ~ ^, and Drr, ~ ^ x {-^f. 
Some additional comments. Our formalism here is not applicable to the non-magnetic 
systems, where M(r) = everywhere. Then we need to determine U in other ways. A 
possibility is utilizing the static numerical linear-response calculations; it gives the infor- 
mation of the static {uo = 0) part of directly (easiest spin-polarization mode at each 
site). Then it will be possible to determine U from such informations together with some 
additional assumptions. In the case of systems like Gd where the d shell and / shell can 
polarize separately, we may need to extend our formulation so as to include non-locality of 
U (e.g. U can be parametrized as Uijki where i,j, k, I are atomic eigenfunction basis for d or 
/ channel). 
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III. RESULT AND DISCUSSION 



A. MnO and NiO 



Fig. [T] shows the calculated SW energies a;(q) for MnO and NiO. (We used 1728 fc-points 
in the BZ for all calculations, including MnAs.) ci;(q) calculated from the LDA is too large, 
as earlier workers have found--. The detailed shape of uj(q) is different from earlier work 

n 

however: in Ref. 'S', peaks in co'(q) occur near 200 meV for NiO, much lower than what we 
find. QSGW predicts uj{q) in good agreement with experimental data. 

The difference of results between QSGW and LDA is understood by Eq. flTB]) . Jrr' 
between nearest AF sites essentially determine the SW energies (exactly speaking, three J 
parameters as shown in Table [T]). The LDA severely underestimates Zi-Eex- This can be 
corrected by LDA+?7, however, Solovyev and Terakurai^ showed that it fails to reproduce 
SW energies as we mentioned in the introduction. This means that the transfer t is also 
wrong in LDA+f/; in fact t is through the hybridization with Oxygen 2p (superexchange). 
In other words, the agreements with SW experiments in QSGW indicates that both of them 
are well described by QSGW. Together with the fact that QSGW showed good agreements 
with optical experimenta^ii^ for MnO and NiO, we claim that our one-particle picture given 
by QSGW captures the essence of the physics for these systems. Our claim here is opposite 
to Refsl5,25 where they claimed that the one-particle picture can not capture the essence. 



B. a-MnAs 

Because a-MnAs is observed to be a FM with a moment of 3.4/1^^, we construct 
assuming a FM ground state. Inspection of the density of states (DOS) in Fig. HJ shows 
that QSGW predicts AE^x ~ 1.0 eV larger the LDA. This difference is reflected in the 
spin moment: Ma= 3.51^b in QSGW, 3.02/iB in LDA. Fig. [2] shows the imaginary part 
of Tr[x+~(q, cu)] along T — K line. Sharp SW peaks are seen at small q; they broaden 
with increasing q. Fig. [3] shows the peak positions, corresponding to SW energies uj{q). 
Hatchmarks indicate the full-width at half-maximum, extracted from data such as that 
depicted in Fig. [2J This corresponds to the inverse lifetime of a SW which decays into spin- 
flip excitations. (Our calculation gives no width for MnO and NiO, because of the large gap 
for the spin- flip excitations). SW peaks are well identified all the way to the BZ boundary. 
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FIG. 1: (color online) Spin wave dispersion a;(q) for MnO and NiO calculated from the LDA and 
QSGW. Solid line without symbols in MnO or dots in NiO (red) are experimental valuea^i^. We 
used experimental lattice constants 4.55 and 4.17 A for MnO and NiO respectively. 
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FIG. 2: (color online) Im part of Tr[x^ (k,a;)] for QSGW. Data are for 6 fc-points, all along the 
F— K line. The fc-point is i/6 K (thus i=6 falls at K). Peak positions and full-width-half-maxima 
are shown in the F— K line of Fig. [3l 

We find that the collinear FM ground state is not stable in the LDA: as Fig. [3] shows, 
a;(q) < around K. (Among all possible collinear configurations, the FM state may be the 
most stable. We did not succeed in finding any collinear configuration more stable than the 
FM one. A similar conclusion was drawn for the PBE GGA functionaF^.) On the other 
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FIG. 3: (color online) Spin wave dispersion w(q) in a-MnAs. QSGW results (circles) are enveloped 
by hatched regions, which indicate the full-width at half-maximum of the spin wave, and is a 
measure of the rate of SW decay. LDA (squares) predicts negative SW energies around K; indicating 
that the collinear FM ground state is not stable. The experimental lattice constants a=3.70 A and 
c/a =1.54 were used. 
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FIG. 4: (color online) quasiparticle DOS for a-MnAs. Lighter hatchings indicate total DOS; 
darker(black and red) hatchings indicate the partial d contribution, whose centers of gravity are 
shown by arrows. The Fermi energy is at zero. 
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hand, QSGW predicts stable collinear ground state, that is, uj{q) > everywhere. However, 
even in QSGW, the SW energies are still low around K, which is a vector that connects 
nearest neighbor Mn sites in x-y plane. If this SW energy is further lowered for some reason, 
we may have a frustrated spin system because of the triangle (honeycomb) lattice of the Mn 
sites. This could be related to the anomalous phase diagram of MnAs, which can easily occur 
through the small changes in lattice structure associated with higher-temperature phases. 

We can qualitatively understand the difference of SW energies between QSGW and LDA 
from the difference of AE^.^. Let us consider the energy difference of FM and AFM states for 



two-site model as illustrated in Ref 
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J. Then the energy gain of a FM pair is independent 
of AEcx when some of majority states are occupied (less than half filling); we measure the 
energy from the majority spin's atomic level as the zero. In contrast, the gain of a AFM pair 
increases with decreasing AEq^. Overall, the LDA with its smaller AEq^, should contain a 
stronger AFM tendency. 



C. Determined Parameters and related Quantities 

Table [I shows the effective interaction (interaction between unit spins). In NiO and 
MnO, f/° as calculated by LDA is much smaller than the QSGW result. This is because 
the LDA underestimates bandgaps in NiO and MnO, thus overestimates the screening. f/° 
is twice larger in NiO than in MnO. This is because Ma{r) is more localized in NiO; in fact, 
the QSGW dielectric constants eoo are similar (eoo = 3.8 (MnO) and 4.3(NiO)^), suggesting 
that the screened Coulomb interaction U{r, r') is similar in the two materials. is smaller 
in MnAs than in MnO, because it is a metal. 

For MnO and NiO, we confirmed that Juji/ is non- negligible only for the three nearest- 
neighbors (NN) (Table [I]). Ji_|_ and Ji_ refer to 1st NN, spins parallel and spins antiparallel, 
respectively. J2 refers to 2nd NN— . Ji+ and Ji_ by QSGW are quite different in MnO, 
while in LDA Ji+ ^ Ji_, resulting in uj{M) ^ in that case. 

For MnAs in QSGW, the expansion coefficients written as {M~^)aa' = — |a;=o, 
is rather dependent on q; nor is {M~^)aa' oc 6aa'- Off-diagonal contributions of {M~^)aa' 
give ~10 % contribution to SW energies. In addition, its inverse of the diagonal element 
l/{M~^)aa is reduced by ~ 0.5fiB at certain points in the BZ. In this case, mapping to a 
Heisenberg Hamiltonian has less clear physical meaning. 
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TABLE I: Magnetic 


parameters calculated by QSGW and LDA 


(in parenthesis). Muffin-tin 


radii R for cations 


were taken to 


be 


2.48 (MnO), 2.33 (NiO), 


and 2.42 (MnAs) a.u. Ma 


is the spin moment 


within the muffin 


tin. Our approximation 


is equivalent to the assump- 


tion for U as 


U{r,r' 




+ LoU^ + ...)ea(r)ea(r') in Eq. 


(l6|). Then Ua is written as 


= Ld'rJ.d'r'ea 


(r)ea(r')C/(r,r', 




= 0). Exchange parameters Ji+,Ji-,J2 si's shown for 


MnO and NiO 


. Total 


1 spin moments for MnAs are 7.00/XB/cell(QSGVF) and S.SD^B/ceh (LDA). 


Our definition of Ji+, 


, Ji_, J2 follows that of Ref. f], except we distinguish Ji_|_ and Ji_^^. 






MnO 




NiO 


a-MnAs 


C/0(eV) 




2.43 (0.95) 




4.91 (1.64) 


1.08 (0.93) 


Ma iflB) 




4.61 (4.35) 




1.71 (1.21) 


3.51 (3.02) 


Ji+ (meV) 




-2.8 (-14.7) 




-0.77 (0.3) 




Ji- 




-4.8 (-14.7) 




-1.00 (0.3) 




J2 




-4.7 (-20.5) 




-14.7 (-28.3) 




Tn or Te (K) 




111 




275 


510 


(experiment) 




122"^ 




523'^ 


400 



D. Calculation of Tn and Tc based on the Heisenberg model 

From obtained Jrr', we estimated Tn (Tc for MnAs) for QSGW (Table [T]) using the cluster 
variation method adapted to the Heisenberg model^^, which assumes classical dynamics of 
spins under H. In NiO, the calculated Tn is only ~ 50% of experiment. There are two 
important effects that explain the discrepancy: (a) QSGW overestimates the dd exchange 
splittingi^ii^, and (6) the classical treatment of quantum dynamics of spins under Ti. Both 
effects will increase Tn. Considering that QSGW well reproduces SW energies (Fig. [1]), the 
errors connected with (a) would not seem to be so serious in MnO and NiO. (6) can be rather 
important, especially when the local moment is small. This is a general problem as discussed 
in Ref. O: Heisenberg parameters that reproduce SW energies well in NiO do not yield a 
correspondingly good Tn. If we multiply our classical Tn by a factor S{S + 1)/S'^ ~ 1.86 
(as 25 = Ma = 1.71), which is the ratio of quantum to classical Tn in mean field theory, we 
have better agreement with experiment. This is what Hutching et al. used^^. On the other 
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hand, evaluation of the quantum Heisenberg model using a Green's function technique show 
that the mean-field theory rather strongly overestimates quantum corrections^. Also, Tn is 
already close to the experimental value in MnO. This is explained in part because correction 
(6) is less important in MnO, since S is larger. Further, we have large contributions to Tn 
from Ji± in MnO, but not in NiO. Around Tn, Ji+ and Ji_ will tend to approach some 
average value, which reduces ci;(q) and therefore Tn (recall uj(M)=0 when Ji+=Ji_). The 
temperature-dependence of J is not accounted for here. 

JfiRi exhibits long-ranged, oscillatory behavior in MnAs: its envelope falls off as |R — R'|^ 
as predicted by RKKY theory for a metal. Consequently, it is not so meaningful to estimate 
Tc from just a few NN, as was done recently^i^^. Shells up to 25th-neighbors are required 
to converge to within 5% or so. The calculated is llOK too high in comparison 
with experiment. Taking (a) into account will improve the agreement; however, there are 
many factors that make a precise calculation very difficult. We also need to take (6) into 
account; in addition, other factors such as assumptions within the Heisenberg model, may 
give non-negligible contributions. 

In conclusion, we present a simple method to calculate spin susceptibility, and applied 
it in the QSGW method. SW energies for MnO and NiO are in good agreement with 
experiments; in a-MnAs the FM ground state is stable, which also agrees with experiment 
(to our knowledge, no SW energies have been published in a-MnAs). LDA results come 
out very differently in each material. By mapping to the Heisenberg model, we estimated 
Tn or Tc. We found some disagreement with experiments, and discussed some possible 
explanations. 
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APPENDIX A: STATIC J{q) CALCULATION HEISENBERG MODEL 

We derive the linear response to an external magnetic field B for the Heisenberg model, 
whose Hamiltonian is given as 

'^^ = — X! X! -^TaT'a'Sxa " Sx'a' + QfJ'B ^ Sxa " Bxa, (Al) 
Ta X'a' Ta 

where Sxa is the spin at Ta (T is for primitive cell, a specify magnetic site in a cell). 
JTaTa = 0. JxaT'a' = -^T'a'Ta- The equation of motion —ihS^a = [H, Sxa] is written as 

^Sxa = Sxa ^ I 2 X! ^TaT'a'^T'a' — Of^B^Ta 1 (A2) 
\ T'a' ) 

We introduce gfis^ = 2b, and Sxa = + ^^xa- i^ the static spin configuration. 
Then Eq. f lA2p reduces to 



^^Sxa — Sl}.^ ^ I 2 X! ^TaT'a'^Sx'a' 1 + ^Sxa ^ I 2 ^ JTaT'a'^T'a' 1 " 2SJ],^ X bxa 
\ T'a' / \ T'a' / 

= ^ (28^^ JxaT'a') X '^Sx'a' ^ ( 2 ^ ^XaT'tt'S^'a' ] ^ ^Sxa " 28^^^ X bxa (A3) 
T'a' \ T'a' ) 

Introducing the Fourier transform, Z\Sxa = ;^ Z^k ^Sa(k)e*^*^'^+'^\ Eq. flASP reduces to 

^^S„(k) = Y. (2S°J,,,(k) - (^2X: 4a"(0)S°„j X ZiS,,(k) - 2S° x b,(k).(A4) 

Assuming ZiStt(k) have 

E + S° J„.Kk) - (^X: •/aa"(0)S°„) X ZiS.,(k) = S° X b,(k). (A5) 

Let us consider the coUinear ground state. Then S° = S^e^ [Sa is the size of spin, including 
sign). We have 

E f^ll^) ^Sa'(k) + E (-^-'(k) - (E ;^'/aa"(0)5a") ^aa') 6, X ZiS,,(k) = 6, X b,,(l<)\6) 

a' \ / a' \ \ a" / / 

Using S = 5+^^^ + 5-^2;^ + 5^e^, and e, x (e^ ± iey) = ^z(e^ ± ie^) we have, 

E(^-^-'(k))5j;(k) = 6+(k). (A7) 

E + Jaa'{\^)] ^a^(k) = ^^(k), (A8) 



2Sa 
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where 



Jaa'(k) = ^aa'(k) " Jaa" (0) Sa'^ 5aa' (A9) 

Only the difference between Jaa'iX) and Jaa'iX) are diagonal parts. These are determined 
so that / d^kJaaCk) = 0. Eq. (|A7| is the same as Eq. (ITill . 
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